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The problem of ferromagnetism - associated with a ground state with maximal total spin - is 
discussed in the framework of a hard-core model, which forbids the occupancy at each site with 
more than one particle. It is shown that the emergence of ferromagnetism on finite square lattices 
crucially depends on the statistics of the particles. Fermions (electrons) lead to the well-known 
instabilities for finite hole densities, whereas for bosons (with spin) ferromagnetism appears to be 
stable for all hole densities. 
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I. INTRODUCTION 

It is generally believed that ferromagnetism should be 
explained in the framework of a Hamiltonian for strongly 
coupled electrons. Such a Hamiltonian on a lattice has 
been proposed 40 years ago by Hubbard. ^-^-^ It contains 
a nearest neighbour hopping term Ht and an on-site 
Coulomb repulsion (of strength U) Hjj for the electrons. 

Both terms are "blind" to the spin of the electrons and 
it is therefore far from being obvious that such a model 
should lead to ferromagnetism. The limit U —^ oo - here 
called "hard-core condition" - enforces that each site can 
be occupied by at most one electron. In this limit, the 
number Q of electrons which can be accomodated on a 
lattice of N sites is limited Q < N. In the charge sector 
Q = N every spin configuration of the electrons is a pos- 
sible ground state. On a square lattice this degeneracy is 
lifted in the sector Q — N — 1 with one hole. Nagaoka^ 
constructed the ground state and found that the total 
spin S{Q) of the electrons is maximal: 

S{Q) = |, Q = N-i (1.1) 

a clear signature for ferromagnetism. 

There has been an intensive search for ferromagnetism 
in the t — J model with iV^ = N — Q holes. Numerical 
results on finite clustersSii at vanishing exchange cou- 
pling {a = 0, see Ea. l|2.1|l ) show that the ground state 
does not have maximal spin except for the Nagaoka case 
Q = N —1. The instability of the Nagaoka ferromagnetic 
state in sectors with 2 and more holes has been investi- 
gated with variational ground states of the Gutzwiller 
^ypgj7j8j9jio .^g^g found that there exists a critical hole 
density 6 = N^/N = 1 — p such that the ground states 
cannot be ferromagnetic for S > 5 crit.- Improvements of 
the variational ansatz led to a successive decrease of 5 crit. 
for the 2D case 

5cr,t. = 0.49i, 0.41iS, 0.292iia, 0.25lii. 

On the other hand Barbieri, Riera and Youngi^ have ar- 
gued that the ferromagnetic ground state emerges in the 



thermodynamical limit N oo, if the hole density van- 
ishes (5^0). The existence of a lower bound 5 > 5 crit. 
for the hole density 5, where ferromagnetism is not pos- 
sible, is indeed a quite general feature, which has been 
found also in higher dimensions and for various types 
of latticesi^ An introduction to ferromagnetism in the 
Hubbard model can be found in Ref . [l 4| . 

In this situation one might ask for the reason why hard- 
core models with fermions (i.e. electrons) fail to "ex- 
plain" ferromagnetism on finite clusters (except for the 
1 hole case). We want to demonstrate in this paper that 
this failure can be traced back to the anticommutation 
relations for the fermion operators. For this purpose, we 
substitute in the t — J model the fermionic degrees of 
freedom by bosonic ones. In order to facilitate the com- 
parison of the fermionic and bosonic version, we assume 
that the bosons carry here as well spin 1 /2 and experience 
the same antiferromagnetic interaction: 

Therefore the Hamiltonians for a hard-core model with 
bosons (6) and fermions (/) are both blind with re- 
spect to the spin. Nevertheless, we find pronounced 
differences (and similarities) in the ground state ener- 
gies Ei{Q, Si{Qj), i — f,b and ground state spins Si{Q), 
i = f,b. 

(a) For Q = 0, 1, iV — 1 the ground state energies and 
ground state spins coincide 

Et{Q,St{Qj) = EfiQ,SfiQ)) (1.2) 
Sb{Q) = Sf{Q)^Q/2 (1.3) 

This means in particular that the construction of the 
Nagaoka ferromagnet state holds for fermions and for 
bosons. The restriction to the one hole sector — 1, 
Q = — 1 is crucial. Moreover one observes a hole- 
particle symmetry for Q = 1, 

E,{Q,S,iQ)) = E,{N^Q,S,iN-Q)) i = f,b. 

(1.4) 

(b) For Q = 2 the ground state energies coincide [ l|1.2|) 
for Q — 2], but the ground state spins are different 

SfiQ^2) = Sb{Q^2)^l, (1.5) 
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i.e. in the bosonic case the ground state spin is maximal. 

(c) For Q = 3, . . . , iV — 2 ground state energies and 
ground state spins are different for the fermionic and 
bosonic version. In the bosonic version the ground state 
spin Sb{Q) ~ Q/2 is maximal as shown in Appendix 1X1 
The ground state energies Eb{Q,Sb{Q)) show the hole- 
particle symmetry (|1.4|) for these Q-values. This symme- 
try is not present in the fermionic version. In the latter 
case the ground state spin Sf{Q) is "erratic" 

We think that these properties are interesting enough 
to justify a detailed study of the ground state properties 
of the bosonic t — J model and to compare its proper- 
ties with the corresponding properties in the fermionic 
version. 

The outline of the paper is the following: 
In Sectiortnlwe review the definition of the t — J model 
and point out the differences between the fermionic and 
bosonic version. Consequences for the ground state en- 
ergies on the smallest (2 x 2) cluster are discussed in Sec- 
tior lllll The phase diagram of the bosonic t~ J model in 
dimensions D = 2 and D = 1 is discussed in Sectior llVI 
andlVl respectively. 

II. THE BOSONIC VERSION OF THE t J 
MODEL IN TWO DIMENSIONS 

Let us first recall the definition of the t—J Hamiltonian 
with fermionic degrees of freedom: 

H = t[Ht+aHj^ a^J/t. (2.1) 

The hopping term 

^* = E E i.^t{^)cAy) + h.c.) V (2.2) 

(x,y) <y 

and the spin exchange part 

H.J = (s(a;)S(y) - }-^n{x)n{y))v (2.3) 

can be expressed in terms of creation [c^{x)) and anni- 
hilation {ca{x)) operators for the electrons, which obey 
anticommutation relations 

{c+(x),c+(x')} = {c,{x),c,,{x'))={) 

(2.4) 

n„{x) = c^{x)cc{x) (2.5) 

is the number operator for an electron at site x with 
spin a. Owing to the "hard-core condition" which is im- 
posed by the projector V double occupancy is forbidden 
i.e. X](T^o'('^) ~ ^T^- The latter can be derived from a 
Hubbard modeli with infinite on-site Coulomb repulsion 
{U/t^oo). 



The construction of a state with Q electrons 

= ctM)---4,A^Qm (2-6) 

by application of creation operators to the vacuum |0) de- 
mands an ordering of all sites on the 2 dimensional square 
lattice. Owing to the anticommutation rules H2.4I) a dif- 
ferent ordering (e.g. with x'[^ , x'^ interchanged) leads to 
a state which might differ from the former one in sign. In 
the following we will denote the traditional model (|2.1|l - 
(|2.3|) with fermions as "fermionic t—J model" . 

Let us now turn to the "bosonic t — J model" , which we 
simply define by substituting the anticommuting creation 
and annihilation operators 

c+(x) a+(x) Ccr{x) -> a„{x) (2.7) 

by commuting ones: 

[a+(x),a+(x')] = {a„(x),a„,(x')\ -0 
[a+ (x) , a„f (x')\ = bcr.a' 5x,x' ■ 

(2.8) 

The eigenvalues of the number operator 

n„{x) = a^{x)aa{x) (2.9) 

are restricted again to nc,{x) =0,1 due to the "hard- 
core condition". This model has been discussed first in 
references 

Quite recently a generalized Hubbard model with 
fermionic and bosonic degrees of freedom has been pro- 
posed and investigated in order to study a mixture 
of ultracold bosonic and fermionic atoms in an optical 
lattice.^ 

The construction of a state with Q bosons by appli- 
cation of creation operators on the vacuum, however, 
is symmetric under the permutation of sites - due to 
the commutation relations H2.8|) . There does not exist a 
"sign problem" . This is the only difference between the 
fermionic and bosonic version of the t—J model [cf. Ref. 
(16)]. One can easily verify, that the action of the spin 
exchange term (|2.3|l is the same in the fermionic and 
bosonic version. On the other hand the hopping term 
Ht acts indeed in a different way on the fermionic and 
bosonic states. In the latter case Ht can be expressed in 
terms of nearest neighbour permutations P(x, y) which 
interchange a particle and a hole at sites x and y: 

Ht ^ P{x,y){n{x)nh{y) + n{y)nh{x)j . 

(2.10) 

Here n(x) — n^(x)+n^{x) = 0, 1 is the number of bosons 
at site X. 

In Appendix ^ the Hamiltonian (|2.10|) is proven to 
have a ferromagnetic ground state for all Q. Moreover 
it is shown that the ground state energies of the bosonic 
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t — J model at a = are symmetric under the particle- 
hole transformation 

Eb{Q,a^Q,N) = EbiN-Q,a^O,N) (2.11) 

as can be seen in Fig. ^ 
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These states are simultaneous eigenstates of the hopping 
part Ht and the spin coupling part Hj 

tHt\F,Q) = EFiQ)\F,Q) (2.16) 
Hj\F,Q) = (2.17) 

such that the eigenvalues Ep{Q) for the t — J Hamilto- 
nian (|2.1|) are independent of a. This is a consequence 
of the fact that the hopping part Ht is blind with re- 
spect to the spin of the electrons. There is no difference 
in the hopping of spin-up and spin-down particles. On 
the other hand, the ferromagnetic eigenstates \F, Q) are 
totally symmetric under any permutation of spin-up and 
spin-down particles. Therefore, the application of the 
nearest neighbour spin exchange coupling 



hj[x,y) = {P[x,y)-l)n[x)n{y) (2.18) 



yields zero: 
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hj{x,y)\F,Q) - 0. 



(2.19) 



FIG. 1: Ground state energies per site of a periodic 4x4 
and a vTO x vTO lattice for both fermionic and bosonic t — J 
models at a = 0. 

Here we also show the ground state energy for hard- 
core fermions, where the hole-particle symmetry ()2.11f) 
is lost (except for the Nagaoka case Q = — 1) 



III. GROUND STATE PROPERTIES OF THE 
FERMIONIC AND BOSONIC t J MODEL ON A 
2x2 PLAQUETTE 

On a plaquette with 4 sites (2 x 2) the t— J Hamilto- 
nian can be diagonalized analyticallyt^ It is interesting 
to compare the ground state energies (Q) and total 
spins S''^P^((5) in the sectors with Q spin-1/2 particles for 
the fermionic and bosonic version. It turns out 



Ef{Q, a = 0, iV) ^ Ef{N -Q,a^Q,N). (2.12) (1) they are the same in the sectors Q = 0, 1, 3, 4 



Moreover the ground state energies of the hard-core 
bosons are below those of the hard-core fermions 

Eb{Q,a^Q,N) < Ef{Q,a^Q,N) .. (2.13) 



Comparing the results from lattices L x L, L = v 10, 4 we 
find that the finite-size dependence is rather small and 
smooth in the bosonic case, in contrast to the fermionic 
case. The result presented in Fig. ^were obtained with 
a Lanczos algorithm. Concerning the fermionic version, 
they agree with Ref. (s*) (Riera, Young). 

In both versions (fermionic and bosonic) the t — J 
Hamiltonian H2.ll) conserves the total charge 



'3 = X! ('^+(^ 



and total spin 



(2.14) 



(2.15) 



such that the eigenvalues of — S{S+\) can be used to 
characterize the eigenstates of the Hamiltonian. Ferro- 
magnetic eigenstates Q) have maximal spin S = Q/2. 



E^P\{)) = 

Ef{l) . 

£;(p)(3) = 

£;(p)(4) = 




g{p) - 

{F) S^f : 
{B) S'-P^ : 
(C) S'(P) : 

Sip) 





1 
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(3.1) 
(3.2) 



(3.3) 
(3.4) 



3to (a > 0) 

where the a intervals (F), (B), (C) are defined by 

3^ 



(C) 



< a < aF(^) = 0.262 

(!) 



a > 2 



< a<2 



Ground state energies with index "F" correspond 
to ground states with maximal spin S = Q/2. In 
the sector with three particles {Q — 3) one finds the 
same energy in the fermionic and bosonic version 
of the t — J model for all a > 0. 
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(2) they are different in the sector Q = 2 
a) In the fermionic version, 



a>0, S'(P)(2)=0 



(3.5) 



the ground state is nondegenerate and has total 
spin 0. 

b) In the bosonic version: 



4"^ (2) 



£;(p)(2) 



0<a<aF(i) = V2, 4''^(2) = 1 

(3.6) 



{\), S^^\2) 



2 

a > ap 



(3.7) 



we find first a ferromagnetic ground state with 
maximal spin 1 and for larger a values degenerate 
spin ground states. 

Note that the ground state energy of the fermionic 
model (|3.5|l is below that of the bosonic model H3.6|l . 
(|3.7|) for a > 0. The two energies coincide at a = 0. 

It is remarkable to note, that some of these similarities 
can be found on all square lattices L x L, L = 2, \/T0, 4. 
Most important is the equality 

Ef{Q = N-l,a,N) = Eh{Q^N-l,a,N) 

N ^ L^, a > (3.8) 

in the 1 hole sector for all a > 0. The a dependence is 
shown in Fig. |21 

In the sector with 2 particles (Q = 2) at a = the 
ground state energies are the same, but the ground state 
spin is always Sf{Q) = in the fermionic, but Sb{Q) = 1 
in the bosonic version, as already stated in (|1.5|l . 



IV. THE PHASE DIAGRAM OF THE BOSONIC 
t - J MODEL IN TWO DIMENSIONS 

In this paper we are mainly concerned with the ques- 
tion, whether the ground state with Q (spin-1/2) parti- 
cles has maximal total spin S = Q/2. The answer to this 
question depends on the strength a of the spin exchange 
part (|2.1|l which prefers (for a > 0) antiferromagnetic 
ordering with total spin S — 0. 

Owing to the property (|2.17|) the expectation value of 
the spin exchange part 
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{E{Q,a)\Hj\E{Q,a)) (4.1) 



FIG. 2: a-dependence of ground state energies in the 1 hole 
sector of a L X L lattice for a > 0. Bosonic and fermionic 
i — J ground state energies for L = vTO, 4 coincide. 



can be considered as an order parameter; it vanishes if 
the ground state \E{Q^ a)) in the sector with Q particles 
has maximal spin S = (5/2, but is nonzero in all other 
cases. 

We will demonstrate for the bosonic version of the t— J 
model on L X L clusters that the ferromagnetic regime: 

ej(p, a) = for < a < apip ^ Q/N) (4.2) 

is bounded by a curve a — apip — Q/N) depending on 
the hole density S = 1 — p. 

In our numerical calculations of the ground states 
\E{Q, a)), Q = TV - 1, iV - 2, ... in the bosonic t - J 
model with periodic boundary conditions, we only found 
ferromagnetic states \F) = \Ef{Q)) for a small enough. 
Fig. 131 shows the order parameter ej{p,a) on a 4 x 4 
cluster in the charge sectors Q = 3, 4, .., 15. 

The phase boundary apip = Q/16) can be read off 
from those a- values where e,j{p,a) changes first from 
to a nonvanishing value. Obviously, ap{p) is monoton- 
ically decreasing with p and vanishes for p = I. Its p- 
dependence is shown in Fig. 2| 

For increasing values of Q, we observe more and more 
jumps in the order parameter ej{p,a). E.g. for Q = 8 
we find three of these jumps. Each of them signals the 
transition to a new ground state with increasing values of 
a. The ground states differ in their total spin S which is 
denoted by the integer {Q even) and halfinteger {Q odd) 
numbers at the branches between two jumps. E.g. for 
g = 8: 



5 = 4 for 0.0 < a < 1.1 = aF(l/2) \Ep) 

5 = 2 for 1.2 < a < 1.6 \Eb) 

S = 1 for 1.7 < a < 2.4 = ao(l/2) \Ec) 

5 = for a> 2.4 1^;^) 

So far we were not able to localize the spin couplings 
for the ground state with total spin 5 = 3, which we 
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FIG. 3: Order parameter £j(p, a) for a A'^ = 4 x 4 cluster with 
charge sectors Q = 3, 15. The order parameter for each 
Q is turning away from zero ai a = arip = Q/N). Small 
numbers denote the total spin of the bosonic t — J ground 
states. 



expect to emerge in the interval 1.1 < a < 1.2. For 
increasing values of p — Q/N , ap{p) approaches zero and 
it is more and more difficult to localize the a values for 
those ground states with higher spin S — Q/2—l,Q/2 — 
2, ... We performed higher resoluted searches with step 
width Aa = 0.001 and found e.g. in the Q = 15 sector 
the ground state with total spin = 1.5 at a = 0.077. 
In contrast, a similar search for the S* = 2.5 ground state 
in the Q = 7 sector failed. 

As an alternative we also tried to localize those a val- 
ues: 



ao(p), P 



Q 
Q 



even, S — 
1 



P= Q odd, 2 

where the total spin S of the ground state is first minimal. 
The curves ao{p), ai (p), and apip) are shown in Fig. 01 
All curves are monotonically decreasing with p. We 
have looked for the finite-size dependence of the curve 
ap (p) by comparing results on L x L lattices for 
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FIG. 4: p-dependence of the lines ao{p), Qi/2(p), otp^p) for a 
iV = 4 X 4 cluster. 



On the larger systems with L > 4 we assumed that the 
ground state spin Sh = Q/2 is maximal and the curve 
ap{p) was extracted from those a- values, where the spin 
is lowered first by one or two units. As can be seen from 
Fig. [SI the finite-size effects can be rather well accounted 
for by an effective charge density 
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such that ap{p') scales with the system size N = L x L. 
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FIG. 5: Finite-size effects of arip) by comparison of results 
on L X L lattices for L = VW, 4, ^18, \/20, \/26. The finite- 
size effects has been accounted for hy p' — p — j^. 

A similar finite-size analysis for the curves ao(p) and 
ai/2{p), which mark the boundaries of the antiferromag- 
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netic phase with minimal ground state spin demands sim- 
ulations on much larger systems. They might be acces- 
sible with Monte Carlo techniques if it is possible to de- 
termine the ground state spin with high accuracy. The 
question of physical interest is, wether the transition re- 
gion between the ferro- and antiferromagnetic phase in 
Fig. ^shrinks to zero in the thermodynamical limit. 

V. COMPARISON OF THE BOSONIC AND 
FERMIONIC t J MODEL IN ONE DIMENSION 

In one dimension the difference between the bosonic 
and fermionic version of the t~ J model can be absorbed 
in the boundary conditions. It turns out that the ground 
state energies with periodic boundary conditions on a 
= 16 site system coincide for all a if Q is odd 

Eb{Q,S,iQ),a,N) = Ef{Q,S,iQ),a,N). (5.1) 

For Q even one observes slight differences for a ~ 2, as 
can be seen from Fig. |H| 




a 



FIG. 6: Comparison of the bosonic and fermionic ground state 
energies per site {Ebj{N,Q,a)/N) for a one dimensional pe- 
riodic chain A'^ = 16 (D = 1) for charges Q = 2,3, 4, 5, 6. 

At a = the equality (|5.1(l strictly holds for all Q- 
values and all the ground state energies satisfy the hole- 
particle symmetry (|1.4|l . 

Concerning the total spin of the ground state Si{Q) 
i — b, f at a = we observe 

a) a maximal value Sb{Q) = Q/2 for all Q 

b) a maximal value Sf{Q) — Q/2 only for odd Q 

For Q even the ground states appear to be degenerate 
with different values for Sf{Q) of the total spin. This 
leads to the "erratic" behaviour observed already in the 
two dimensional system D — 2. 

Differences between the bosonic t — J model for D = 1 
and D = 2 become apparent in the variation of the order 



parameter (|4.1|l and the total spin Sb{Q) with a. They 
are shown in Fig. [7|for D = 1 and in Fig. Ofor D ~ 2. 




0.5 1 1.5 2 2.5 
a 



FIG. 7: Order parameter ej{p,a) for a A*' = 1 x 16 cluster 
with charge sectors Q = 3, 15. Small numbers denote the 
total spin of the bosonic t — J ground states. 

The main difference between the two cases D = 2 and 
D = 1 \s the following: 

The ferromagnetic domain < a < ai?(p) - where the 
total spin is maximal - shrinks to zero in the ID case. 
In particular for Q > 3, the boundary ap{p) turns out 
to be below 0.1 already on finite chains. 

VI. DISCUSSION AND PERSPECTIVES 

In this paper we have demonstrated that the emergence 
of ferromagnetism in a hard-core model crucially depends 
on the statistics of the constituent particles, i.e. whether 
they are fermions or bosons. We studied the t — J model 
in the (traditional) fermionic (/) and in a bosonic (6) ver- 
sion and compared ground state energies Ei{Q, Si{Q), a) 
and total spins Si{Q,a), i — f, b, in both models. In par- 
ticular, we found for a = that - in contrast to hard-core 
fermions - hard-core bosons (with spin 1/2) have a ferro- 
magnetic ground state with maximal spin Sb{Q) ~ Q/2 
in all charge sectors Q. 

This means that the corresponding Hamiltonian Ht 
(|2.10|) leads intrinsically to a ferromagnetic interaction 
between spin- up and spin-down particles, as we have 
demonstrated first on finite square lattices L x L, L ~ 
2, y/TO, 4 numerically. A general proof is given in Ap- 
pendix ^ 

The ground state energies for hard-core bosons and 
fermions behave differently on finite lattices. They show 
hole-particle symmetry (I2.11|) in the bosonic version, 
whereas this symmetry is not present in the fermionic 
version. We do not know whether this breaking is a finite- 
size effect and whether the ground state energies (Fig. ^ 
converge to each other in the thermodynamical limit. If 
we switch on the spin exchange coupling a [ (|2.3|l . (|2.18|) ]. 
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the bosonic t — J model in two dimensions (D = 2) devel- 
ops a ferromagnetic regime Q < a < apip), where ap{p) 
is monotonically decreasing with p. For D — 1, apip) 
shrinks to zero. 

Ferromagnetism is observed in the fcrmionic and 
bosonic t — J model in the sectors Q = 1 (1 particle) and 
Q = iV— 1 (1 hole). Here, the two t~J models are unitary 
equivalent, which implies that the sign problem for the 
electrons can be absorbed in a unitary transformationn^^ 
In other words "electrons behave as bosons". This be- 
haviour can be found as well for hard-core electrons in 
ID with periodic boundary conditions and Q odd. 

On one hand all ground state energies Ef{Q,a — 0) 
Q = 1, .., N — 1 show the hole-particle symmetry ()1.4|l 
and are strictly identical with the corresponding quan- 
tities Ef,{Q,a — 0) in the bosonic version. The ground 
state spins Sf{Q) are maximal for Q odd, but this is not 
the case for Q even, where the ground state appears to 
be degenerate with different total spins. 

Let us finally comment on possible more realistic hard- 
core boson models for ferromagnetism. On one hand 
these bosons have to carry spin, on the other hand the 
spin statistics theorem demands integer spin for bosons. 
It is quite straightforward to extend the hopping Hamil- 
tonian l|2.1()|l to spin 1 bosons. In this case, the number 
of particles at site x: 

n{x) = ni(x) + n_i(x) -I- no(a^) = 0, 1 (6.1) 

has to be identified with the sum of spin 1,-1,0 particles. 
Such a Hamiltonian is again invariant under particle- 
hole transformations. It defines a hard-core boson model 
with 4 degrees of freedom at each site: 3 for the spin-1 
particles, 1 for the holes. 

APPENDIX A: THE GROUND ST ATE O F THE 
HOPPING HAMILTONIAN (ITTnil 

We want to present first a general proof that the fer- 
romagnetic state with all spins up (A^+ = Q) is indeed a 
ground state of Ht (|2.10l) in the sector with Q particles. 
Here, the permutation operator P{x, y) for a spin (-I-) 
and a hole can be represented: 

P{x,y) = i(l + r(x)f(y)) (Al) 

in terms of Pauli matrices at neighbouring sites x and y 
which act on particle x+ ^-nd hole Xh states as 

n{x)xj{x) ^ {2n+{x) - l)xj{x) , j^+,h 











, Xh = 





It is easy to verify that (|2.1U|) in the sector — Q can 
be mapped on the XX spin model in two dimensions: 

^* = (^i(^)^i(j^)+^2(a:)r2(2;)). (A3) 



The Hamiltonian IjASp is invariant under the inversion of 
all "r spins" (HU 

UHtU+ = Ht U = Y[n{x) (A4) 

X 

which explains the hole-particle symmetry H2.11|l . 

Note, that all the matrix elements in Ht have a nega- 
tive sign. Therefore we can apply the Perron-Frobenius 
theorem)^ which states that the ground state iPf{Sz = 
Q/2) is nondegenerate and that in the defining basis all 
components are greater than zero. 

Since Ht conserves the total spin H2.15|l . the applica- 
tion of lowering operators 

S- = ('^i(^)-*'^2(y)) (A5) 

onto the ground state ippiSz — Q/2) with eigenvalue 
Ep{Q) yields eigenstates of Ht : 

(5_)>f (5. = Q/2) = (V;(Q/2)|5^^!!|^(Q/2))-i/2 X 
V'f(^. -Q/2-n) 

(A6) 

Hti'FiS, = Q/2 -n) = EF{Q)ijF{S, = Q/2 - n) 

(A7) 

in the sector with 5"^ = Q/2 — n with the same eigenvalue 
Ef{Q)- Indeed these states arc again ground states in the 
sector 5"^ = Q/2 — rt for the following two reasons: 

i) The state iPf{Sz = Q/2 — n) has again only positive 
components in the defining basis in the sector Sz = 
Q/2-n. 

ii) The matrix elements of Ht in the sector Sz = Q/2— 
n are all negative, such that the Perron-Frobenius 
theorem can be applied again: The ground state is 
nondegenerate and has only positive components. 
Therefore iPf{Sz = Q/2 — n) must be the ground 
state for S'^ = Q/2 — n. 

Finally, we want to stress that all the states iPf{Sz — 
Q/2 — n) are symmetric under all permutations: 

(^P{x, y) - l) {n+ {x)n^ (y) + n_ {x)n+ {y)j x 

i^F{Sz=Q/2-n) = 

(A8) 

for the quantum numbers of the particles at sites x and 
y, which guarantees that 

S^t^FiSz^Q/2-Ji) - SiS+l)^F{Sz^Q/2^n) 

S^Q/2. 

(A9) 
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